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ABSTRACT 

We present a 3D dynamical model of the orbital induced curvature of the wind-wind 
collision region in binary star systems. Momentum balance equations are used to 
determine the position and shape of the contact discontinuity between the stars, while 
further downstream the gas is assumed to behave ballistically. An archimedean spiral 
structure is formed by the motion of the stars, with clear resemblance to high resolution 
images of the so-called "pinwheel nebulae" . A key advantage of this approach over grid 
or smoothed particle hydrodynamic models is its significantly reduced computational 
cost, while it also allows the study of the structure obtained in an eccentric orbit. The 
model is relevant to symbiotic systems and 7-ray binaries, as well as systems with 
0-type and Wolf-Rayet stars. 

As an example application, wc simulate the X-ray emission from hypothetical 
0+0 and WR-f O star binaries, and describe a method of ray tracing through the 3D 
spiral structure to account for absorption by the circumstellar material in the system. 
Such calculations may be easily adapted to study observations at wavelengths ranging 
from the radio to 7-ray. 

Key words: hydrodynamics - mcthods:numcrical - stars:early-type - X-rays:stars - 
starsibinarics - stars:winds 



1 INTRODUCTION 

Colliding winds occur in various types of stellar binaries, in- 
cluding those with massive OB and Wolf-Rayet (WR) stars, 
, lower mass eruptive symbiotic systems containing a white 
dwarf and red giant star which undergo a "slow nova" out- 
burst, and binary systems which contain one or two pulsars 
blowing a pulsar wind(s). 

High-spatial-resolution observations are revealing many 
interesting features in such systems. In massive O-f O and 
WR-I-O binaries, radio interferometry has spatially-resolved 
emission from non-t hermal electrons at the apex of the wind- 
wind collisi o n (e.g., Williams et al.l Il997l : iDoughertv et al] 
I2OOOI . I2OO5I : IContreras et all 120011 ). Beautiful "pinwheel" 
struc tures which trace dust emission can also be observed 
(e.g.. iTuthill et al.lll999l . I2OO6I . I2OO8I : iMonnier et al] 1 19991 : 



iMarchenko et al.ll2002i r The shape of these structures can 
be described by archimedean spirals which are believed to 
follow the wind-wind collision region in systems where the 
winds are of very unequal momentum. 

Colliding winds also play a key role in eruptive symbi- 
otic systems, where a hot, fast, diffuse wind from a white 
dwarf companion interacts with a slow massive wind from 
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a Mira type primary staiQ. The class of eruptive symbiotics 
can be divided into two further subtypes: classical symbi- 
otics, in which the bolometric luminosity remains constant 
and outbursts typically last about 100 days (Z And is an 
example), and the more powerful eruptions known as sym- 
biotic novae, where the bolometric luminosity increases by 
a factor of order 10-100 on a timescale of about a year, and 
the system stays in an active state for > 10 yrs (well-known 
examples are V1016 Cyg, HM Sge, and AG Pegasi). 

Colliding winds may also play a key role in the 
newly discovered class of systems called 7 -ray binaries 
j Aharonian et all l2005al lbl: [Albert et aLll2006l ). The nature 
of these systems is still controversial, though in the case 
of PSR B1259-63, it is clear that a relativistic wind from 
a pulsar collides with the stellar wind from a Be star. 



^ Such systems are to be distinguished from the interact- 
ing wind phenomenon which occurs in AGB binaries where 
density structures in the AGB wind are created either 
due to the reflex-action of t he evolved star around the 
centre-of-mass of the system |Mastrodemos &: Morrid Il999t 
iMauron fc Hugginsl[2006l: [H^|2007^. or due to gravit ational fo- 
cusing }Gawrvszc2ak~MikoIaiewska &: Rozvczkal I2OO2I ) . Here we 
focus exclusively on systems involving the interaction of winds 
from separate stars. 
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The orbit is highly eccentric (e = 0.87), and has a period 
of 3.4 yr (jjohnston et al.ll2005f ). The TeV 7-ray emission 
arises from the Inverse Compton cooling of uhra-relativistic 
elect rons accelerated at the pulsar wind termination shock 
fe.g. iKhangulyan et al.|[2007l ). In contrast, the nature of the 
sources LS5039 and LS I +61 303 is less clear, since the 
type of compact object ha s not been established beyond 
doubt (iRomero et all l2007l: iDubus. Cerutti fc Henril [20oi : 



IKhangulyan. Aharonian fc Bosch-Ramonll2008l ). 

While there has been much progress in mod- 
elling the dynamical structure of the colliding winds 
in early-type binary systems, the majority of work has 
been limited to 2D (e.g 



19921: iGayley Qwocki fc Cranmerl l 19971 : 
Pittardlll998l l2007l : IZhc kov 2003). 3-di 



Stevens. Blondin fc PoUockl 



Pittard et al]|l997l : 



nttardlllblijai . ll^UUYI : IZ;nc kov 2mji). a-dimen sional hydrody - 
namical calculations ha ve be en performed bylPittardI (| 19991 ). 



IWalder fc Folinil (|2002l ). and lLemaster et al.l (120071'). w hile a 
ballistic model was •Resent ed bv l Harries iniTm). An 
SPH model has recently been computed by lokazaki et al.l 
l|2008l ). Dyn amical models for symbi otic novae have been 
presented bylcirard fc WillsonI |l983) and lKenny fc Tavloj 
(|2005l . I2OO7I ), while 3D hydro d ynam ical models have been 
presented by IWalder fc Folinil (2OOOI ). Models of the wind- 



wi nd collision in classical symbiotics have been presented 
by iMitsumoto et all (|2005l) and J3isikalo ct al. ( 20_oi). Rel- 
ativistic hydrodynam ics (|Bogovalov et al] boosl ) and SPH 
l|Romero et al.1 120071 ) models have been used to investigate 
the wind-wind collision in pulsar wind binary systems. 

Although dramatic improvements in computational 
power and techniques in recent years have spurred the de- 
velopment of 3D models of colliding winds, such work re- 
mains computationally expensive, and it is still difficult to 
perform simulations of CWB's even on high performance 
parallel machines when the orbital eccentricity is high. We 
therefore present a new method which captures the flow dy- 
namics while requiring less computational resources. 

At its heart, our approach adopts the equations for the 
ram pressure balance between the two winds as detailed by 
ICanto et al.l (| 19961 ). In this work it is assumed that both 
winds are highly radiative, rapidly cool, and fully mix. While 
these assumptions are only relevant in close binaries, it pro- 
vides a convenient starting point and the position of the 
contact discontinuity is unlikely to drastically change even 
if the wind-wind collision is essentially adiabatic. Then, at 
some distance downstream of the apex of the wind-wind 
collision region (WCR) the flow is assumed to reach a ter- 
minal speed and to thereafter flow ballistically (i.e. no net 
force acting upon it). This ballisti c treatment has similar- 
ities to many previous works (e.g.. 



iGirard fc Wi llson 19871: 

iHarries et aLll2004l: [Kenny fc Tayloill2007l ). The derivation 



Canto et al. 1 19961 ) has also been widely used to model 



obser vable properties (e.g.. lFoellmi et al ] |2008l : lHenle7 et al.l 
I2OO8I ). 

This paper is organised as follows. In § [2] we explain 
the steps necessary to construct our dynamical model of the 
wind- wind collision. §|3] shows how it can be used to simulate 
the X-ray emission and circumstellar absorption arising from 
the WCR in early-type binary systems, though this is but 
one example of the potential use of such a model. In § |3] we 
summarize and conclude our findings, and outline possible 
future directions. 
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Figure 1. Schematic diagram showing the location of the shocks 
bounding the CD for varying values of the cooling parameter. 
Left: both winds have X ^ 1 s-nd ^^e therefore largely adiabatic. 
Right: The primary wind has x ^ 1 s-nd is strongly radiative, 
whereas the secondary's wind is adiabatic. The scenario on the left 
represents the O+O and WR-l-0 binaries investigated in this pa- 
per, whereas the scenario on the right represents what is thought 
to occur in 77 Car . 



2 THE DYNAMICAL MODEL 
2.1 Overview 

In the model the orbit is calculated in the frame of one of the 
stars (herafter referred to as the primary star). The winds 
are assumed to reach their terminal speeds before they col- 
lide. The contact discontinuity (CD) is split into two sec- 
tions, to account for the effect of orbital motion: 

i) A region close to the apex of the WCR where the 
flow from the stagnation point is accelerating along the CD 
(hereafter caUed the "shock cap"). The shock cap is termi- 
nated where the flow is assumed to become ballistic (the 
"ballistic point"), the exact point being calibrated against 
hydrodynamical models (see 12. 3|) . While the properties of 
the shock cap are assumed to be axisymmetric, orbital mo- 
tion introduces an aberration angle which means that the 
symmetry axis and the line of centres of the stars are not 
colinear (see ^ I2.2|l . 

ii) A region beyond the ballistic point where the flow 
along the contact discontinuity is unaffected by the primary 
and secondary stars' gravity, ram pressure from the winds, 
or thermal pressure in the WCR. If the stellar winds have 
differing speeds, the flow in this region is assumed to move 
with the speed of the slower wind, since this is the wind 
which responds least to the orbital motion of the stars, and 
dominates the absorption in the system (in the models pre- 
sented in this paper, both winds have the same speed of 
2000 kms~^ - see Table [5}. This region is termed the "bal- 
listic CD". 

By separating the CD into these two sections we can 
model the effect of the winding of the CD around the stars 
and the subsequent absorption by the un-shocked winds. We 
do not attempt to model the shocks which bound either side 
of the CD in this work, as in many circumstances the shocked 
gas efficiently cools and is compressed by the ram pressure 
of the pre-shock wind into a thin dense sheet coincident with 
the CD. For instance, in symbiotic novae, the hot wind is 
likely to be strongly radiative (see Fig. 4 in I Kenny fc Tayloil 
[iooi), as of course is the cool wind, and our model therefore 
gives the position of the shocked gas and the dense spiral 
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Figure 2. Ratio of the cooling to characteristic flow timescale for 
the hot shocked gas in a coUiding winds binary as a function of 
orbital period. Cooling becomes important once X ^ 3- 



shells which subsequently form. Strong radiative cooling is 
also a feature of the WCR in many massive binaries. In 
7] Car , for example, the primary LBV wind is so dense (and 
slow) that it is str ongly radiative around the entire orbit 
(|Pittard et aLlllOgSl ). The importance of cooli ng in the WCR 
can b e quantified using the cooling parameter jStevens et al.l 



^cool 
tesc 



M_7 ' 



(1) 



rfi 



where vs is the wind velocity in units of 1000 km s" 
is the separation of the stars in units of lO^^cm, M7 is the 
mass-loss rate of the star in units of 10"^ M0 yr~^ , fcooi is 
the cooling time, and tosc (~ d/cs, Cs is the postshock sound 
speed) is the characteristic time for hot gas near the apex of 
the WCR to flow downstream. In practice, hydrodynamical 
simulations show that the wind collision region (WCR) is 
adiabatic for x ~ 3, whereas for x ~ 3 it cools rapidly 
(Fig.©. 

Fig. [2] shows the value of x as a function of orbital pe- 
riod for each shocked wind in a hypothetical O-l-O star bi- 
nary with a circular orbit, wind speeds of 2000 kms~^ , and 
mass-loss rates of 10~® M0 yr~^ and 2 x 10~^ MQyr~^ for 
the primary and secondary star of masses 50 and 30 Mq 
respectively. Clearly both of the shocked winds are largely 
adiabatic, even down to an orbital period of 10 days (in 
shorter period systems the stars are close enough together 
that acceleration/deceleration of the winds needs to be con- 
sidered). However, in a hypothetical WR-fO system where 
the mass-loss rates of the primary and secondary stars are 
now 2 X 10~^ Mq yr"^ and 10"** Mq yr"^ , and both stars 
have masses of 50 Mq and wind speeds of 2000 kms^^, cool- 
ing is important for orbital periods < 1 yr. If the WR star is 
a WC subtype, cooling is important for periods up to several 
years, sin ce cooling is more e fficient with such abundances 
(see, e.g.. IStevens et aLlll992l ). Thus, Fig. [2]shows that the 
denser winds from WR stars are likely to produce radiative 
shocks in many instances, though the O-l-O systems will usu- 
ally be adiabatic unless the orbital period, P < 10 d, or the 
winds are slower and/or denser than assumed above. 

The postshock winds of both the primary and secondary 




Primary star 



Figure 3. Schematic diagram describing the wind-wind interac- 
tion between the two stars represented by circles. 



stars in the simulations discussed in § [3] are largely adia- 
batic. In such cases, the temperature of the hot gas in the 
WCR as a function of distance downstream from the stagna- 
tion point at the apex of the WCR has been determined by 
iKennv fc Tavlon (|2005l ). With this information it is possible 
to derive the width of the post-shock layer, and hence the 
position of the shocks, as a function of downstream distance. 
However, this is beyond the scope of the present work. 

In the following sections we detail the modelling of the 
shock cap and ballistic CD. 



2.2 The shock cap 

The shape of the shock cap is determined from momentum 
balance requirements. The surface density and velocity of 
the fl ow along the shock cap are obtained from Eqs. 29 and 
30 of ICanto et al.l (| 19961 ) (the latter scaled to the speed of 
the slower wind). Assuming the winds are already at their 
terminal velocity when they reach the shocks, the locus of 
the CD is R{9ci), where R is the distance from the centre 
of the primary star and 6ci is the angle between the vector 
to the primary star and the line-of-centres (see Fig. [3| . The 
ratio of the wind momenta is given by 



M2V0 



(2) 



where Mi, Vooi, M2, and Voo2 are the mass- loss rates and 
terminal velocities of the primary and secondary stars re- 
spectively. 

The shock cap is symmetrical about the line of centres 
before the effects of orbital motion are introduced. The 2D 
(r, z) coordinates of points on the shock cap in units of the 
stellar separation, dscp, are 

tan^ci 



tan^c2 + tan^ci 
z tan 6c2 ■ 



(3) 
(4) 



To determine the coordinates in 3D, the 2D arms of the 
WCR can be rotated azimuthally. The x,y, and z vectors 



cap , Vcap , 



and 



respectively) from the center of the 



primary star to coordinates on the shock cap are then 



Xcap = dscp {z cos Lu — r sirnu cos (^), 
ycap = ciscp(2:sina; + rcostJCOsC), 
Zcap = dscp(rsm(), 
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Figure 4. Plan view of tlie shoclc cap at various orbital phases 
in a system with an orbital eccentricity, e = 0.9. The calculations 
are performed in the frame of the primary star (located at the 
origin and not to scale) with the secondary star orbiting in an 
anti-clockwise direction. As the separation of the stars increases 
the linear size of the shock cap also increases in direct propor- 
tion. This is most noticeable when comparing the shock cap at 
(p = 0.0 (periastron, red) and <j) = 0.5 (apastron, turquoise). The 
aberration (skew) of the shock cap due to orbital motion is not 
considered in this plot. 



where uj is the true anomaly of the orbit and C, is the az- 
imuthal angle subtended between a coordinate on the sur- 
face of the shock cap, the line of centres, and the xy (orbital) 
plane. 

The number of coordinate points on the shock cap is de- 
termined by the values of Scioo, S9ci and S(. With 56ci = 1° 
and = 18° (i.e. 20 azimuthal points per 2D rz value), 
the shock cap consists of ~ 10'^ separate coordinate points. 
Eq. 29 of lCanto et al] (| 19961 ) is used to determine the tan- 
gential velocity along the CD, and thus the position of the 
ballistic point in 2D axisymmetry. 

Since the size of the wind-wind collision scales with the 
orbital separation, dramatic variations occur in systems with 
highly eccentric orbits, as shown in Fig. |4] where e = 0.9; the 
high eccentricity means that the shock cap at periastron has 
a linear scale which is 20 times smaller than that at apastron. 
Such high eccentricities occur in two of the most well-known 
colliding winds systems, r; Car and WR 140, and also in PSR 
B1259-63, one of the 7-ray binaries. 

Another effect resulting from orbital motion is the aber- 
ration (skew) of the apex of the WCR due to the net velocity 
vector of the orbit (i.e. the motion of the secondary star rel- 
ative to the primary star). The skew angle, /x, which is the 
angle between the symmetry axis of the shock cap and the 
line of centres of the stars is approximated by 



tan /i : 



(5) 



where the speed of the slower wind is used. In the frame of 
the primary star, 

1/2 

(6) 



«orb = 



G{Mi + Ma) ( - i 

' Ob. 



Figure 5. Variation of the skew angle of the shock cap due to 
orbital motion for eccentricities of e = 0.3 and e = 0.9, with 
P = 1 yr and Mi + M2 = 80 Mq . The angle peaks at phase 
= 0.0, when the stars are at closest approach and their relative 
orbital velocity is a maximum. 



The aberration is significant in symbiotic novae because 
of the low wind speed of the cool star (for instance, a sym- 
biotic system with e = 0.0, Mi + M2 = 2.5 M©, and 
rfscp = 10 an has i^orb = 15 kms^^ , which is comparable 
to the speed of the cool wind). In contrast, the aberration 
is small in early-type binaries (for instance, an O-l-O binary 



with e = 0.0, Ml -I- M2 = 80 Mq, and dsop = 4.3 au has 
VoA = 130 kms^^ , which is typically much smaller than 
the wind speeds), unless the orbit has high eccentricity. In 
such cases the magnitude of the skew varies throughout the 
orbit. 

Fig. [S] shows how [i, varies throughout the orbit for the 
Model A O-l-O star binary with parameters as in Table [2] 
and with e = 0.3 and 0.9. A peak value is reached at peri- 
astron passage ((/> = 0.0) when the relative orbital speed of 
the stars reaches it's highest value, and the lowest value of 
occurs when the stars are at apastron (0 = 0.5) and the 
relative orbital velocity is a minimum. The variation of [s, 
between apastron and periastron increases with the eccen- 
tricity of the orbit. The skew angle can affect the proximity 
of regions of the shock cap to the primary star around pe- 
riastron (Fig. [Sjl , and the resulting level of occultation and 
attenuation. 



2.3 The ballistic CD 

To construct the large-scale 3D structure of the WCR, gas 
packets are released from the endpoints of the shock cap at 
specific phase intervals with a velocity equal to the slower 
wind, Wai. The x,y, and z components of the velocity of gas 
leaving the end of the shock cap at a specific orbital phase 
are given by: 



for stars of mass Mi and Afa and an orbital semi-major axis, 
a. 



Vx ~ Val COs(tJ — fl + \ COS C, ) , 

Vy — fsi sin(a; — /i -|- AcosC), 
Vz = Vsi sin 
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Figure 6. Plan view of tiie sliock cap at = 0.0 for e = 0.9, 
P = 1 yr and M1 + M2 = 80 Mq, witli (green) and witliout (red) 
taking account of orbit induced skew. The secondary star orbits 
in the anti-clockwise direction. Note that the shocks either side 
of the CD are not displayed. 



d = VB.(tc-%) d = VA.(t(,-tA) 




Figure 7. Schematic diagram showing how position vectors to 
points on the ballistic CD are obtained. If the secondary star is 
currently at position C, points pi on the ballistic CD are found 
by determining the position and flow direction of gas leaving the 
end of the shock cap when the star is at position B and advecting 
the flow by the time difference between these two orbital phases. 
This process is repeated until the position of the ballistic CD has 
been traced back the desired number of orbits (in this work the 
structure is traced back over 2 orbits). Distances of lOO's-lOOO's 
of au can easily be covered, as required. 
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Figure 8. 3D plots viewed from above the orbital plane showing 
the position of the contact surface separating the stellar winds 
from a hydrodynamic simulation (top) and from the dynamical 
model with a varying cut-off percentage for the transition between 
the shock cap and the ballistic CD (bottom). The shading in 
the hydrodynamic simulation identifies the different stellar winds, 
where a value of 0.5 marks the CD. In these simulations an orbital 
period of 19.7 days, eccentricity, e = 0.0, stellar masses. Mi = 
AI2 = 30 Mq, terminal wind speeds of fool = foo2 = 1500 kms-i 
were used. The mass-loss rates of the primary star was Mi = 
3 X lO~^M0yr~^ and that of the secondary star was Mi = 
1 X IO-'^Mq yr-i . 



where A — OcIod is the asymptotic half-opening angle of the 
contact discontinuity viewed from the star with the stronger 
wind. 

The baUistic part of the CD is then constructed by 
considering a sequence of previous positions of the ballistic 
points at the termination of the shock cap, and the current 
position of the gas flow from these points given that they 
move along linear trajectories (see Fig. (Tjl . 

The position of points on the ballistic part of the CD 
(a;cD, J/CD, and zcd) at the time t is given by their position 
at the time they were emitted from the end of the shock 



cap (xcapil/cap, and Zcap) plus the distance they have since 
travelled at velocity Vx,Vy,Vz, i.e. 

a;cD (t) = a;,ap {t~T)+ {t - T)T, 
ycD (t) = j/cap {t-T)+ i,y {t - T)T, 

ZCD (t) = 2cap {t~T)+ {t - T)T, 

where T is the time elapsed since the flow left the end of the 
shock cap. 

The number of coordinates in the ballistic CD is de- 
pendant on the number of phase steps around the orbit, 
the number of orbital revolutions followed, and the number 
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Figure 9. 3D plots viewed from above the orbital plane showing the position of the shock cap (red) and the structure of the ballistic CD 
(green) at various orbital phases from a simulation with e = 0.3 (sec Table [2] for the other relevant parameters). The skew to the shock 
cap due to the orbital motion of the secondary star is included. From top left to bottom right: <f> = 0.0, 0.1, 0.5, andO.7. The curvature of 
the ballistic CD at = 0.0 is caused by the high velocity of the secondary star around periastron passage. Note that the shocks either 
side of the CD are not displayed. The white area downstream of the ballistic CD in the corner of each plot is an artifact of the plotting 
software. 



of azimuthal steps (i.e. SQ. In this work, the ballistic CD 
consists of 2000 coordinate positions along each azimuthal 
trajectory (1000 per orbit traced). 

Tests performed using a 3D hydrodynamics code con- 
firm that the Coriolis force, which causes the curvature to 
the WCR, becomes significant once the flow from the stag- 
nation point is accelerated to 70% - 90% of the terminal 
speed of the slower wind (Fig. [S| , and the gas is at a dis- 
tance from the stars of order the stellar separation. Both 
of these conditions are satisfied by the 85% cut-off attained 
via calibration of the dynamic model against hydrodynamic 
models. Interestingly, varying the value of the cut-off per- 
centage has the effect of improving the fit to one spiral arm 
but reducing the quahty of the fit to the other arm. Using 
the Model A paramters (Table[2]) the off-axis distance of the 
ballistic point from the line of centres, rmax, increases by a 
factor of ~ 3 between 70 % and 90 % (Table [l|, whereas the 
opening angle of the shock increases by roughly a half with 
a more linear relation. In 13.31 we show that there is little 



Table 1. Transition points between the shock cap and the bal- 
listic CD for varying percentages of the postshock flow velocity, 
''max is calculated using Eq|4] fcloo is the opening angle of the 
shock at the transition point measured from the primary star (see 
Fig.©. 



Cutoff 


^max 




(%) 


(dscp) 


(°) 


70 


0.53 


33 


80 


0.78 


40 


85 


0.99 


44 


90 


1.40 


49 



difference in the X-ray lightcurves when this percentage is 
varied slightly. 

Fig- m shows the effect of the motion of the stars on the 
ballistic CD on scales of the order of the semi-major axis. 
The curvature of the CD close to the end of the shock cap 
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is greatest when the relative orbital velocity of the stars is 
high. The smooth connection of the ballistic CD to the shock 
cap indicates that the assumptions inherent in the model are 
good at this level. 

The structure of the ballistic CD at large scales is shown 
in Fig. [To] for a range of orbital eccentricities. At low or- 
bital eccentricities, the s piral structures re s emble the 30 hy- 
drodynamical models of IWalder fc FoUnil (|2000l , l2002h and 
iLemaster. Stone fc Gardineil ( 200?! ) , t he dust spira l mod - 
els of the pinwheel neb ula WR 104 by iHarries et al] (j2004l ) 
and lTuthill et al.1 (|2008|). and the CW o model for symbiotics 
developed bv iKennv fc Tavlo j (|2007h . Note, however, that 
this figure shows the projection of the CD onto the orbital 
plane, and not the position of the shocks either side of it. If 
the shocked region were largely adiabatic, the shocks would 
stand off from the CD and the width of the spiral structure 
on the orbital plane would be somewhat greater. 

At e = 0.9 the secondary star moves very quickly 
through periastron, resulting in the projected CD (which 
encompasses the region of unshocked secondary wind) thin- 
ning to the left of the stars. In contrast, there exists a large 
region of unshocked secondary wind to the right of the stars, 
as the secondary star moves slowly around apastron. This 
creates a low density cavity in the primary wind. The X- 
ray attenuation in such systems will depend on the orbital 
phase, as well as being sensitive to the position of the ob- 
server, and in principle may vary widely. For instance, in a 
system like Car, the primary wind is very dense and much 
more strongly absorbing than the secondary wind. An ob- 
server located on the positive x-axis at infinity will predom- 
inantly view through the low density unshocked wind of the 
secondary star, whereas an observer on the negative i-axis 
will predominantly view through the high density unshocked 
wind of the primary star. As the column density scales di- 
rectly with the density of the gas, these observers will see sig- 
nificantly different X-ray lightcurves. On the other hand, if 
the primary wind is more rarefied than the secondary wind, 
this behaviour reverses. 

Finally, we note that in systems with highly eccentric 
orbits, the amount of attenuation at phases around apas- 
tron may depend on the skew angle of the shock cap which 
occurs around periastron. This is because the skew angle 
of the shock cap affects the duration and phase where pri- 
mary/secondary wind material is emitted in a certain direc- 
tion. Depending on the viewing angle into the system, the 
inclusion of aberration effects may result in a variation in 
the attenuation to emission concentrated near the apex of 
the shock cap due to the alteration in path length through 
the more strongly absorbing wind. 



3 AN EXAMPLE APPLICATION - X-RAY 
EMISSION AND ABSORPTION IN AN 
EARLY- TYPE BINARY 

As an example application of the model described in § [2l 
we consider the X-ray emission from hypothetical 0-|-0 and 
WR-fO-star colliding wind binaries. Model A is an 0-1-0 
binary with an orbital period of 1 yr and semi-major axis 
a = 4.3 au. Model B examines the increasing effects of ab- 
sorption as the orbital period is reduced to 1 month. The 
third model system (Model C) consists of a WR star with 



a mass-loss rate of 2 x 10~^ Mq yr~^ . The high velocities of 
the stellar winds are sufficient to cause the postshock gas to 
emit at X-ray wavelengths, and both winds are essentially 
adiabatic in the systems considered (see Table [2]) . For the 
three models considered we use a distance of 1 kpc, ISM 
column of 5 X 10^^ cm^^, and orbital eccentricity, e = 0.3. 



3.1 The X-ray emission 

The X-ray emission from the WCR is a function of 
the gas temperature and density. Since the dynamical 
model discussed in the previous section does not con- 
tain such information, we use a grid-based, 2D hy- 
drodynamical calculation of an axis-symmetric WCR to 
obtain this. The numerical code is second-ord e r ac- 
curate in time and space (|Falle fc Komissarovl Il996l : 
iFalle. Komissarov fc Joardej Il998t ). The resulting emission 
as a function of off-axis distance is then mapped onto the co- 
ordinate positions in the 3D dynamical model. In this way 
we obtain the benefit of effectively modelling the thermo- 
dynamic and hydrodynamic behaviour responsible for the 
production of the X-ray emission, while simultaneously ac- 
counting for the effect of the motion of the stars on the 
large-scale structure of the WCR and the subsequent wind 
attenuation. Since the hydrodynamic calculation is 2D, the 
computational requirements remain low. 

The X-ray emission calculated from each hydrodynamic 
cell in the WCR is r{E) = n^VA{E,T), where n is the 
gas number density (cm~"^), V is the cell volume (cm'^), 
and A{E, T) is the emissivity as a function of energy E 
and temperature T for optically thin gas in coUisional 
ionization equilibrium (erg cm'^s^^). A{E,T) is obtained 
from look-up tables calculated from the M EKAL plasma 
code iLeidahl. Osterheld fc Goldsteinlll995l , and references 
therein) containing 200 logarithmically spaced energy bins 
in the range 0.1-10.0 keV, and 101 logarithmically spaced 
temperatures from 10* to 10^ K. Solar abundances are as- 
sumed for the O-star winds and the WR wind is assumed to 
have WN8 abundances (mass fractions of: H/He=0, C/He 
= 1.7 X 10"*, N/He = 5 X 10"^ and O/He = 1 x 10""). The 
emissivity of solar abundance gas is shown in Fig.[TT]and the 
corresponding opacity is shown in Fig. [T2] The WN8 emis- 
sivities are very similar to those at solar abundance. Opacity 
values are also similar for solar and WN8 abundances, with 
the most significant difference being a factor of 2 increase at 
10* K at energies below ~ 1 keV. 

The emission values are then appropriately scaled for 
the changing stellar se paration around the orbit (Lx oc d~p 
in the adiabatic limit. fstevens et al.lll99j ) and placed onto 
the 3D shock cap and ballistic CD. Emission values are as- 
signed to points within 3 dsep of the apex of the WCR. This 
accounts for ~ 90 per cent of the 0.1-10 keV emission and 
> 99 per cent of the 2-10 keV emission. 



3.2 The attenuation 

To compute X-ray lightcurves, the orientation of the ob- 
server relative to the system must be specified. Since the 
model assumes the orbit of the stars is in the xy plane, 
viewing angles into the system can be described by the in- 
clination angle that the line-of-sight makes with the z axis. 
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Figure 10. Plots of the large-scale structure of the ballistic CD in the orbital plane at = 0.0 as a function of eccentricity, e = 
0.0, 0.3, 0.6, and 0.9 from top left to bottom right, with the other parameters given in Table[2] The semi-major axis of the orbit is 4.30 au, 
and the orbital period is 1 yr. The black region on the plots shows the projection of the 3D contact discontinuity onto the orbital plane, 
which encompasses the unshocked secondary wind. The white region mostly traces the unshocked primary wind. Note that the shocks 
either side of the CD are not displayed. The variation of the velocity and separation of the stars around the orbit increases with the 
orbital eccentricity, which results in the increasing asymmetry of the projected CD. The tapering of the black region at the tail of the 
spiral (most noticeable in the e = 0.0 and e = 0.3 plots) is due to the fact that the ballistic part is only traced back for two orbits and 
that the leading and trailing arms of the CD finish in different directions. Further details of the interaction at e = 0.3 arc shown in Fig. |9] 
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Table 2. Summary of the wind and orbital parameters of the model systems, r] is the wind momentum ratio, and a the semi-major axis 
of the orbit. In each model the wind speeds and stellar radii adopted were Voci = i'oo2 = 2000 kms~^ and Rti = i?»2 = 10 R©. Cooling 
parameters are calculated for an orbital separation of deep = a. Th e half opening angles, 6cioo, of the WCR are comparable to the 61° 
and 41° for Models A and C respectively calculated using Eq.(3) of lEichler fc Usovl l|l993l 'l. 



i, and the angle the projected line-of-sight makes with the 
major axis of the orbit, 9. Positive values of 9 correspond to 
projected lines of sight in the prograde direction from the 
positive X axis. The components of the unit vector along the 
line-of-sight, u, are thus 



cos 6/ sm J, 
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Figure 11. Emissivity of the solar abundance gas in the spec- 
tral energy range 0.1-10.0 keV at various temperatures (K) as 
calculated from the MEKAL thermal plasma code. 
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Figure 12. Opacity (gcm^^) in the spectral energy range 0.1- 
10.0 keV at various temperatures (K). 



There are 3 ways in which the intrinsic X-ray emission 
can be attenuated. First, it can be occulted by the stars (this 
effect, of course, is greatest in short period systems). Sec- 
ond, there will be some absorption through the un-shocked 
stellar winds. Finally, there will be attenuation through the 
shocked gas in the WCR. The latter is only important in 
systems where the shocked gas of at least one of the winds 
is strongly radiative, otherwise the gas in the WCR remains 
hot and the attenuation through it is small. But if there is 
significant cooling, as for example occurs when the cool wind 
in symbiotic systems is shocked, a thin, dense, and cold layer 
of gas is formed at the CD, which will be a significant source 
of attenuation in the system. Significant cooling of the WCR 
can also occur in early-type binary systems of which r] Car is 
an example. We now describe how attenuation by each of 
the above-mentioned methods is calculated in our model. 



3.2.1 OccuHation by the stars 

An important line-of-sight effect in binary star systems is 
occultation, particularly in the case of eclipsing binaries. To 
calculate this effect in our model, a line-of-sight is traced 



Figure 13. 3D plots showing occultation effects for two different 
viewing angles at periastron {<j) = 0.0) with e = 0.3 and an orbital 
period P = 1 month. The occulted emission from the shock cap 
is represented in green (this is essentially the "shadow" of the 
occulting star - the secondary star is hidden behind the shock 
cap). The top panel shows the case where i = 90° and 8 = 0°. 
The primary star is directly in front of the secondary star and 
the apex of the WCR is occulted. In the bottom panel i = 70° 
and 9 = 0°. The line-of-sight peers over the top of the primary 
star and the "shadow" moves below the apex of the WCR - less 
emission is occulted as a result. The finite number of coordinates 
on the shock cap, together with the trapezium-shaped emission 
regions, account for the particular shape of the "shadow" in this 
case. Aberration of the WCR is included in both plots. 



from each emitting region on the shock cap and ballistic 
CD, and its distance of closest approach to the centre of 
each star is calculated. If this distance is less than the radius 
of the star, and the star is in front of the emitting region, 
then occultation occurs, and none of the emission from the 
emitting region being considered reaches the observer. 

A visual representation of the occultation of the WCR 
by the primary star is shown in Fig. 1131 The degree of oc- 
cultation can be reduced by reducing the inclination angle i 
(since the strongest X-ray emission occurs at the apex of the 
WCR). Although not shown, the phase at which the maxi- 
mum occultation occurs can be altered by changing the value 
oie. 
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Figure 14. 0.1-10.0 keV lightcurves for Model B with e = 0.0, 0.3, 0.6, and 0.9 with occultation included. Left: The variability shown over 
an entire orbital period. Right: The variation over the orbital phase range 0.9-1.1. Occultation is highest at periastron for the e = 0.9 
lightcurve because the size of the shock cap relative to the primary star is smallest at this time. Orbital motion induced skew causes 
the pre-minimum intrinsic luminosity to be lower than the post-minimum, this effect increases with orbital eccentricity. For all curves 
i = 90° and 9 = 0°. Interstellar absorption has not been considered in these calculations. 



Occultation causes little change to the observed lumi- 
nosity over the entire orbit for the Model A system. This is 
due to the relatively small size of the stars in comparison 
to the extended emitting region for the 0.1-10.0 keV X-rays. 
Occultation effects become more noticeable in shorter period 
systems, and/or those with highly eccentric orbits (since the 
linear size of the shock cap is oc 1/dscp). Occultation is also 
favoured where one (or both) of the stars has a large stellar 
radius (e.g. r] Car, Parkin et al., in preparation), and when i 
is large. For instance, the e = 0.9 lightcurve in Fie:. ll4l shows 
a pronounced occultation effect at orbital phase <j) ~ 1.00, 
during which the emission falls sharply by a factor of 2. 
The width of the minimum due to occultation effects is very 
narrow as the high eccentricity means that the stars move 
rapidly through periastron, but the depth of the minimum is 
large (~ 75% of the intrinsic 2-10 keV emission is occulted). 

3.2.2 Absorption by the un-shocked stellar winds 

For inclinations, |i| ^ 7r/2—Scioo, the line-of-sight from emit- 
ting regions near the apex of the WCR will intersect the 
CD numerous times as it spirals out, and thus traverses first 
through one wind and then the other, etc. The total column 
density along a line-of-sight is then the sum of the individ- 
ual column densities along the specific distances travelled 
in each wind. Accurate knowledge of where the line-of-sight 
intersects the CD, and the density of the gas at any point in 
space is therefore required if the total column density along 
a given sight line is to be calculated. 

To determine if and where an intersection through the 
CD occurs, the shock cap and ballistic CD are tesselated 
into a sequence of triangular planar facets constructed be- 
tween three neighbouring coordinates {P_^, P_^, and P_^. To 
determine if the line-of-sight intersects a given triangle the 
normal to the plane in which the triangle lies, n, is calculated 
from 

n = (£b-£J X (P^-PJ. (7) 
The dot product of n with the line-of-sight vector gives the 



angle between the line-of-sight and the plane. If the resultant 
angle is non-zero the line-of-sight vector will intersect the 
plane in which the triangular facet lies at some point in 
space. 

The component vectors to the intersection point (xint, 
yint, and 2int) are found by substituting the line parameter 
at the intersection point, 

^x^cap ~t~ ^^y^cap ~t~ ri^Zcap ^const /qv. 

K = , (o) 

Tlxltx + riyUy + n^u^ 

into the equations 

^int — 3^cap ~t~ f^Ux^ 
Vint = ycap + KUy, 
^int — Zcap ~t~ hiyUz , 

where the equation of the plane with normal n and vector 
components rix, riy, and is 

n^x + Uyy + n^z + riconst — (9) 

In general, the intersection occurs outside of the trian- 
gular facet. Unit vectors are constructed between the corner 
points of the facet and the intersection point to determine 
whether the intersection occurs within its boundaries. The 
three dot products between these three unit vectors gives 
the angles 6a, 6b, and 6c. Only if the intersection point lies 
within the boundaries of the triangular facet will the equa- 
tion 6a + 6b + dc = 2n be satisfied (see Fig. I15|) . By looping 
over the entire sequence of triangles, every possible intersec- 
tion of the line-of-sight with the CD is determined. 

With the coordinates of the intersection points (a;int, 
yint, and Zint) , it is a simple task to calculate the column den- 
sity through the unshocked winds, ctw. Lines-of-sight which 
pass very close to the stars sample the acceleration region 
of the wind. Therefore, we use a /3-velocity law of the form 

/ pa" 

v{r) =Voo[l- ^ j (10) 

to determine the density of the wind at radius r from the 
star. P describes the acceleration of the wind with f3 = 0.8 
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Figure 15. The method used to determine whether a given line- 
of-sight intersects the CD through a specific triangular tessel. 
Triangles are constructed over the entire shock cap and ballistic 
CD from adjacent coordinates and the unit vectors between the 
corners of the triangle and the intersection coordinate are used to 
determine if the intersection occurs within the boundary of the 
triangle. 



appropriate for O star winds l|Lamers fc Cassinellilll999t ). 
Because the width of the WCR is not considered in our 
model, the volume of unshocked wind and the resulting at- 
tenuation are overestimated, though this approximation will 
not have a signifcant impact on our results. 

3.2.3 Absorption by the shocked stellar winds 

As already mentioned, the attenuation of X-rays through the 
shocked wind(s) needs to be considered if one or both winds 
strongly cool. In the O-l-O and WR-|-0-star binaries con- 
sidered in this section, the shocked gas is largely adiabatic. 
However, for completeness, we discuss here a methodology 
for calculating the absorption due to X-rays intersecting a 
cold dense layer of postshock gas at the CD. This is ap- 
plied to models of rj Car in Parkin et al. (in preparation) . In 
Figs. [Tg [m [23 EH [23 US and [21] this effect does not need 
to be considered. 

The surface density, as, of the postshock gas along 
the CD, when both winds have x 1; has been com- 
puted bvlGirard fc Willsoiil (|l987l ). ICanto et al.1 (|l996D . and 
iKennv fc Tavloij (|2005h ." In each of these works, turbulence 
in the postshock flow is assumed to fully mix the material 
from both winds and the surface density calculated is for 
shocked gas from both winds. Alternatively, if only one of 
the winds is radiative (i.e. the other remains largely adia- 
batic), or the postshock flow is assumed not to mix, then the 
surface density can be calculated from consid ering conser- 
vation of mass flux (e.g.. lAntokhin et al.ll200"3 ). To calculate 
th e surface densities i n Figs 1 161 and 1171 we have used Eq.(30) 
of ICanto et al] (| 19961 ). 

In our model, the ballistic part of the WCR is asym- 
metric due to orbital motion. Since the pre-shock flow is 
practically tangential to the CD at this point, we calculate 
the total surface density of the postshock winds (which in 
this subsection are assumed to cool) in this region by con- 
sidering conservation of mass flux. The surface density of 



the postshock gas close to the apex of the WCR varies by 
over an order of magnitude between periastron and apastron 
when e = 0.9. 

Since the width of the cool dense layer of gas along- 
side the CD is not infinitely thin, the degree of absorption 
through it depends on the angle subtended between the line- 
of-sight and the normal to its (i.e. the CD's) surface, 7. The 
column density intersected by the line-of-sight is therefore 

(11) 



cos 7 



where is the actual surface density of the cooled layer. 
When the line-of-sight becomes closely tangential to the CD, 
cTg can become large, even if itself is not particularly large. 

The maximum value of is constrained by the cur- 
vature of the WCR and the finite path length through the 
shocked gas. To determine the maximum path length re- 
quires knowledge of the width of the cooled layer and its 
radius of curvature at the point of interest on the CD. On 
the shock cap the density of the cooled postshock region, 
Pps, can be determined by equating the ram pressure of the 

J ireshock gas with the thermal pressure of the postshock gas 
Kashi fc Sokeill2007l ). 



(12) 



where Tps is the temperature of the cooled postshock gas 
(Tpa is taken to be 10'' K), p(r) is the preshock gas density 
as a function of distance from the respective star, ^ is the 
angle between the preshock velocity vector and the tangent 
to the shock surface, and mn and fee are the mass of a hy- 
drogen nucleus and Boltzmann's constant respectively. The 
thickness of the cooled layer, /shock, is then 

/shock = — — ■ (13) 

Pps 

The thickness of the cooled layer in the ballistic CD region 
cannot be calculated in this manner because the shocks are 
now fully oblique. Therefore, a linear extrapolation is used 
to determine the downstream thickness. 

The radius of curvature at a point on the shock cap is 

ds 

~di ' 

where ds is the distance between two points on the WCR 
and dt is the difference in the unit vectors tangent to the 
WCR at those two points. Consideration of the maximum 
path length through the cool dense layer, dmax, then gives 
the maximum value for — as 



A = 



(14) 



d„ 



Is 



8A 

^shock 



(15) 



The skewing of the shock cap due to orbital motion 
will evoke an asymm etry in the postshock gas density 
(jLemaster et al.l l2007l l . This is naturally accounted for in 
Eas.[T2land[T5] 



3.2.4 The observed emission 

In the hypothetical binary systems considered in this pa- 
per, the shocked gas in the WCR remains hot as it flows 
out of the system and thus contributes insignificantly to the 
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Figure 16. The effect of the various attenuation mechanisms on 
the observed 0.1-10.0 keV emission as a function of orbital phase 
for Model A. The un-attenuated intrinsic emission (I), intrinsic 
emission with occultation (I + O), intrinsic emission with occu- 
latation and wind attenuation (I + O + W), and intrinsic emis- 
sion with occultation, wind attenuation, and attenuation due to 
intersection of the shocked gas (I -|- O + W -|~ SD) are shown. 
Attenuation through the shocked gas (SD) is only important if 
the shocked gas can cool efficiently and is shown here purely for 
illustration of its effect. The assumed viewing angles are i = 90° 
and 9 = 0°, and the other model parameters are listed in Table [2] 
T he skew to the WC R due to orbital motion is included. Eq.(30) 
of lCanto et"al] l ll996l) has been used to calculate surface densities 
as this provides an approximate upper limit to the attenuation by 
shocked gas. Interstellar absorption (~ 5x lO^'^cm"^) is included. 
The variation in the intrinsic emission is proportional to 1/dsop. 



Fig. [16] demonstrates the effect of including the vari- 
ous attenuation mechanisms on the resultant emission. As 
previously mentioned, occultation causes little reduction in 
emission because of the minute size of the stars compared to 
the extended WCR (Fig. I13|) . For the assumed position of 
the observer, absorption by the unshocked winds increases 
as the stars approach each other and reaches a maximum 
at periastron. For illustrative purposes we also show the at- 
tenuation that occurs if the postshock gas cools and forms 
a thin dense layer along the CD (this does not occur in the 
systems considered since the shocked gas remains largely 
adiabatic as it flows out of the system). When the line-of- 
sight becomes closely tangential to the WCR the path length 
of X-rays through the shocked gas and the subsequent atten- 
uation via this mechanism reaches a maximum. This can be 
seen in the small dips in the lightcurve at orbital phases 0.18 
and 0.82, with corresponding peaks in the emission weighted 
column shown in Fig. 1171 The emission weighted column 
density is calculated as Eatot^o/S/o, where the summation 
is over all sightlines to emitting regions and crtot is the total 
column density (cm~^) along each sightUne. This weighting 
is more informa t ive th an the column densities presented in 
iLemaster et al.l (|2007h which were only calculated along a 
single sight-line into the system. 

Resolution tests have determined that the minimum 
number of phase steps required for convergence of the atten- 
uated X-ray lightcurves is dependant on the ratio Uorb/'i^si, 
with of order 1000 phase steps required for an orbit with 
e = 0.9 and Uorb/wsi ~ 1. 



3.3 Results 
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Figure 17. Variation of emission weighted column density as a 
function of orbital phase. Interstellar absorption is not included. 



absorption. Hence the total column density along a given 
line-of-sight is the sum of the column densities through the 
unshocked winds. The attenuation declines as the line-of- 
sight leaves the system, and is negligible at the distances 
which our model extends to (the distance the wind flows 
over two orbits). 

Absorption cross-sections for solar abundance gas at 
10* K are used to obtain the optical depth, r, along specific 
hnes-of-sight in 200 logarithmically spaced bins over the en- 
ergy range 0.1-10.0 keV. The observed attenuated emission, 
lohs = /oe~^, where 7o is the intrinsic emission. 



3.3.1 The X-ray lightcurve 

In this section we compute X-ray lightcurves for the hy- 
pothetical systems considered. For the O-l-O systems we use 
emissivity and opacity data calculated assuming solar abun- 
dances (Figs. [TTl and 1121 respectively) . For the WR wind we 
use data appropriate for WN8 abundances. 

Fig. [18] shows the synthetic lightcurves produced for 
models where the transition between the shock cap and the 
ballistic CD is varied. There is a maximum divergence of 
~ 6% between ceises where the transition occurs at a ve- 
locity cut-off of 70 % and 90 % of the speed of the slower 
wind, which shows that the resulting lightcurves are not very 
sensitive to this assumption. 

Varying the orbital inclination angle changes the 
amount of attenuation that the intrinsic emission suffers 
on its way to the observer. However, there is little circum- 
stellar attenuation for Model A (Fig. 1191 left panel), and 
the synthetic lightcurves are almost identical over the en- 
tire orbital period. This is because the emitting volume is 
large (so occultation by the stars is negligible), and be- 
cause the stellar separation is wide, so that the circum- 
stellar density at the WCR is relatively low. Attenuation 
effects become more prominent if the orbital period is re- 
duced to 1 month (Fig. 1191 right panel) , and distinct differ- 
ences in the lightcurves occur around periastron. The i = 0° 
lightcurve is smooth, and reflects the fact that the increase 
in the intrinsic emission due to the changing orbital separa- 
tion (Lx ex 1/dscp) more than offsets the peak in attenuation 
through the primary wind at periastron. Increasing the in- 
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Figure 19. Synthetic 0.1-10.0 kcV emission lightcurves for various inclination angles, i, for Model A (left, a = 4.3 au) and Model B 
(right, a = 0.81 au). Model parameters are indicated in Table[2] Interstellar absorption is included. 
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Figure 18. Synthetic 0.1-10.0 keV lightcurves for different tran- 
sition positions between the shock cap and the point where the 
flow in the WCR is assumed to behave ballistically for Model A 
(period of 1 year, semi-major axis of 4.3 au). The transition point 
is specified in terms of a percentage of the terminal speed of the 
slower of the two winds, v^i. The assumed line-of-sight is i = 90° 
and 8 = 0°. Interstellar absorption is included. 



clination enhances the attenuation around periastron. The 
dip seen in both the i = 60° and i = 90° curves is offset 
from the time of periastron (0 = 1.0) because of the skew to 
the WCR caused by orbital motion. As aheady mentioned 
in § 13.2.11 the stars fail to provide any significant echpse of 
the emitting region. 

Fig. [20] examines the dependence of the observed emis- 
sion on the angle subtended between the line-of-sight and 
the semi-major axis. The 8 = 90° and —90° curves appear 
to be almost identical reflected copies around orbital phase 
(j) ~ 0.5, with the differences around periastron being due to 
the aberration of the WCR. Absorption does not strongly 
affect the observed emission (even if the orbital period is 
reduced to 1 month) as the density contrast between the O- 
star winds is not very large. The largest difference between 
the model results (~ 25%) occurs at periastron. A compar- 
ison between current observational data and such models 
may allow constraints to be placed on the orientation of 
specific 0+0-star systems. 



I ' ' 1 ' ' 1 1 
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Orbital phase 

Figure 20. Synthetic 0.1-10.0 keV emission lightcurves for Model 
B (see Table |2]l with i = 90° and various line-of-sight angles, 6. 
The observed emission is clearly sensitive to the viewing angle. 
Interstellar absorption is included. 



The higher primary mass-loss rate in the WR+0 sys- 
tem leads to a greater depenence of the observed emission 
on the line-of-sight (Figs. [21] and 122(1 . as well as higher 
X-ray luminosities. The minimums in the curves close to 
periastron in Fig. 1211 especially in the case of the i = 90° 
curve, are the result of the the X-rays passing through the 
dense WR wind. However, there is again little difference be- 
tween the i = 60° and i = 90° curves at apastron as the 
WCR is viewed predominantly through the less dense O- 
star wind, though the attenuation at lower inclinations is 
slightly higher as the apex of the WCR is viewed through 
the denser wind from the WR star. Rotating the line-of-sight 
within the orbital plane again causes significant alterations 
to the observed emission (Fig. I22|) . The 6 = 0° curve sees 
the largest degree of attenuation around periastron and the 
lowest around apastron, with the opposite being true for 
the 9 = 180° curve. As was also the case in Fig. 1201 the 
6 = 90° and —90° lightcurves appear to be almost perfect 
reflected copies of each other. The dip seen in the 9 = 90° 
curve at orbital phase ~ 0.17 marks opposition. Features 
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Figure 21. Synthetic 0.1-10.0 keV emission lightcurves for Model 
C (sec Tablc[2]l with various incUnation angles, i. Interstellar ab- 
sorption is included. 
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Figure 22. Synthetic 0.1-10.0 keV emission lightcurve for Model 
C (see Table [2J with inclination i = 90° and various line-of-sight 
angles, 9. Interstellar absorption is included. 



such as these could be particularly useful for constraining 
the orientation of systems. 

Fig. [23] shows the variation with phase of the emission 
weighted column density for the 3 hypothetical systems con- 
sidered. The column density is highest when viewed through 
the primary wind, and lowest when viewing through the sec- 
ondary wind. It is lowest for Model A, and is approximately 
5x higher when the period is reduced to 1 month (Model 
B). This simply reflects the ~ 5x smaller separation and 
the ~ 30x higher densities. The ~ lOx changes in the col- 
umn density between apastron and periastron in Model A 
curve reflects the 1.86 x change in stellar separation and 5x 
change in wind density (or stellar mass-loss rate) as the line- 
of-sight switches from the secondary wind into the primary 
wind. 

The different slopes of the column density either side of 
periastron are caused by the asymmetry of the WCR. The 
rise in column density at </!) ~ 0.8 begins when the shock cap 
rotates and lines-of-sight start to see the emission through 
the denser primary wind. The rise occurs at an earlier phase 
for the WR-l-0 system because of the lower value of the 
wind momentum ratio and the narrowing of the opening 



1e+20 I ^ ^ ^ ^ ^ 1 1 
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Figure 23. Variation of the emission weighted column density as 
a function of orbital phase for the O-l-O and WR+O-star systems 
considered. The average column is also shown for the WR-I-O- 
star system to demonstrate the difference between this and an 
emission weighted column. Interstellar absorption is included. 



angle of the WCR. At = 0.964 and 0.055 the slope in 
the emission weighted column is reduced, and this feature 
marks the point where the bow shock arms are tangential 
with the line-of-sight. When the emission weighted column 
density is plotted alongside the average column density the 
change in slope occurs at a point where the two curves in- 
tersect. The average column density has a continual rise and 
a peak at periastron. This tells us that the column density 
to the entire emission region reaches a maximum at perias- 
tron, whereas attenuation to the points with highest intrin- 
sic emission remains roughly constant for a short period. It 
is also interesting that both the O-l-O and WR-I-O systems 
with P = 1 year (Model A and C respectively) have flatter 
profiles at maximum column density. This indicates that the 
shape of the column density curve is sensitive to the aber- 
ration and orbital induced curvature of the WCR, and thus 
to the orbital period. 

3.3.2 X-ray spectra 

Fig. [24] shows synthetic spectra at periastron and apastron 
for the simulations performed. The slope of the spectra at 
high energies is the same for the O-l-O systems since the 
preshock velocities, and therefore postshock gas tempera- 
tures, do not change. In all cases spectra at apastron show 
lower flux in the 2.0-10.0 keV energy band, although flux 
below 1 keV is higher. This is because the intrinsic emission 
scales as 1/dsep, but the is weaker when viewed through the 
companion's wind. The low energy turnover in the perias- 
tron spectrum extends to higher energies for the WR-l-0 
system due to the higher mass-loss rate and absorption of 
the WR wind. 



4 CONCLUSIONS 

We have presented a 3D dynamical model of the colliding 
winds in binary systems where both stars drive a significant 
wind. In circular systems, the WCR adopts a spiral shaped 
structure similar to those observed in massive binary star 
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Figure 24. Synthetic 0.1-10.0 kcV spectra at periastron and 
apastron for Model A (top), B (middle), and C (bottom). The 
viewing angles are i = 90° and 8 = 0°. Interstellar absorption 
is included. The corresponding lightcurves for these spectra are 
shown in Figs.[l9l [20] and [21] 



systems. In systems with eccentric orbits, the shape of the 
WCR becomes increasingly deformed as the eccentricity in- 
creases, with the winds increasingly being channeled into 
a specific direction. A major advantage of the model is its 
low computational cost and the fact that it can be easily 
adapted to model a wide variety of observational data (from 
the radio to 7-ray) and systems (from early type binaries, 
to 7-ray binaries with a pulsar wind, to symbiotic novae). 

As an example exercise, the X-ray emission from hypo- 
thetical 0+0 and WR+O-star systems was modelled. The 



intrinsic emission was computed from a 2D grid-based hy- 
drodynamical model of the WCR, and then mapped onto the 
surface separating the winds in the 3D dynamical model. 
Absorption due to the unshocked stellar winds (and also 
cooled postshock material) can be considered, although in 
the hypothetical systems that were modelled the gas in the 
WCR remains largely adiabatic as it flows out of the system 
so that only the former is calculated. Ray-tracing through 
the 3D spiral structures then gives the attenuated emission, 
and synthetic spectra and lightcurves are produced. 

The lightcurves and spectra show that observational 
characteristics of the X-ray emission from early-type bina- 
ries can be reproduced. For instance, the model with a 1 
year orbit (Model A) is representative of w ide 0-1-0 binaries 
such as HD 15558 (|De Becker et all 120061) . and in this par- 
ticular system could be useful in determining whether there 
are two or three counterparts. The results from the 1 month 
orbit simulation (Model B) are instead most applicable to 
X-ray observations of clo se O-l-O binarie s such as l Ori- 
onis (jPittard et al.l [ioooh. CvgOB2#8A (|De Becker et all 
l200d ). and HD 93403 (jRauw et all I2OO0I ). to name but a 



few. The model can also tackle systems with different abun- 
dances for each wind such as WR-l-O-star systems. Our 
WR-I-O star model with a 1 year orbital period (Model 
C) is ap plicable to systems like WR25, WR108, WR133, 
WR 138 (|van der Huchtl200ll ). and among the WN stars and 
WR19, WR125, WR137, WR98a, WR104, a nd WR140 
(|Pollock et al.ll2005l : IPittard fc DoughertvllioO^ among the 
WC stars. 

Mass-loss rate determinations can be made from com- 
parison of t he predicted magnitude of the X-ray flux with ob- 
servations (|Stevensll99^ : |Pittard fc CorcoranH2002l ) . In prin- 
ciple it is possible to use the shape of the X-ray lightcurve 
to constrain the inclination and orientation of the system. 
Our results reveal that for wind momentum ratios of order 
0.2, the lack of significant absorption means that this will be 
very difficult if applied to O-fO-star systems with periods of 
order one year, but becomes possible for orbital periods of 
order one month. The variation in absorption is much more 
significant when the wind momentum ratio is lower and the 
density of the winds is more disparate. This is the case for 
WR-fO, LBY+O, and LBV+WR systems. 

In future work we will apply the dynamical model to 
the X-ray and forbidden line emission from rj Car, the X-ray 
lightcurve of WR 140, and emission line profiles of colliding 
wind binaries. 
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